Ground state fluctuations in rung-dimerized spin ladders 
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Treating an exactly rung-dimerized spin ladder as a reference model we study perturbatively 
zero temperature quantum fluctuations in spin ladders with slightly destroyed rung-dimerization. 
Analytic expressions are obtained for the gas parameter (density of rung-triplets) and the ground 
state energy per rung. At a strong diagonal frustration as well as at a rather strong antiferromagnetic 
rung coupling these results well agree with the previous numerical calculations. 

PACS numbers: 75.10.Pq 



I. INTRODUCTION 

Perturbation theory for rung-dimerized spin ladders 
first considered in seminal worksi*^ was then developed 
in a number of papers^i^i^. In all of them an array of iso- 
lated rungs was treated as the exact zero-order reference 
model. The main virtue of this approach is a possibility 
to perform calculations up to a very high order. From the 
other hand its principal defect is triviality of dispersion 
in excitation spectrum of the reference (isolated rungs) 
model. According to this obstacle it seems reliable that 
obtained results will be well applicable only to narrow 
band spin ladders with AEmagn ^ -E-gap (Egap is a gap, 
AEmagn IS a magnon bandwidth). However this condi- 
tion may be failed for real compounds. For example a 
recent neutron scattering experiment^ in the spin ladder 
compound (C5Di2N)2CuBr4 gives AEmagn/Egap « 0.72 
while for the spin ladder compound La6Ca8Cu2404i this 
ratio is^ about 5.3. 

In this paper we suggest a more general perturbative 
approach to magnetic excitations in spin ladders by treat- 
ing the exactly rung-dimerized spin ladder— iSiiiS as a refer- 
ence model . As in the isolated-rungs casei'^'^'^'^ the cor- 
responding ground state still is an array of rung-dimers 
(singlets) 
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(|0)„ is a singlet state associated with n-th rung), how- 
ever now an excited triplet (denoted by |1)„) can move 
along the ladder from one rung to another. This process 
entails a nontrivial dispersion for an elementary excita- 
tion of the reference model. The latter will be called a 
rungon because in fact it is an excited single rung coher- 
ently propagating inside the rung-dimers bulk. It will be 
implied that the perturbed model belongs to the rung- 
dimerized phase, so its ground state is a dilute gas of 
excited rungon pairs. In this picture an elementary ex- 
citation magnon should be considered as a bare rungon 
interacting with the perturbed ground state. 

The main technical difficulty of the suggesting ap- 
proach is non integrability of the reference model 
Indeed as it will be shown in forward an evaluation 
of the lowest (second) order perturbative correction to 



the ground state energy utilizes only two-rungon wave 
functions^'^" of the reference model. The next, i. e. one- 
magnon level requires a knowledge of the three-rungon 
spectrum. Due to the non integrability this problem is 
of a principal mathematical difficulty. Nevertheless for 5 
integrable cases the three-magnon problem may be solved 
by the Bethe Ansatz^°. 

In this paper the manifested approach is applied for an 
evaluation of two main ground state parameters, namely 
the gas parameter (or equivalently the density of rung 
triplets) 



p= lim -(0|g|0). 
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and the normalized energy density 
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Here H, N and Emng are correspondingly the Hamilto- 
nian, a number of rungs and the excitation energy of an 
isolated rung. The operator Q ~ J2nQn, whose local 
density satisfies the following relations 
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is the number operator for rung-triplets. Probably it is 
better to denote the gas parameter by n (concentration) 
instead of p (density). However in this paper the index 
"n" is utilized especially for enumerating of ladder rungs. 

Both the quantities p and e are of interest for ascer- 
taining of boundaries of the rung-dimerized phas o^^'^^ . 
Also p may be considered as the governing parameter for 
the dilute gas approximatio n^^i^^i^^ . The parameter e 
was usedi^ for evaluation of the spin ladder entropy. 

For a future development of the perturbation theory 
both p and e may be considered as the corresponding 
governing parameters. Indeed it is crucial to obtain a sys- 
tem of conditions on the coupling constants under which 
the perturbative approach is valid. In fact it should be 
a requirement for some relevant dimensionless governing 
parameters to be small. Usually only p was suggested on 
this role. Namely the mean field Green function approach 
to rung-dimerized spin ladders may be developedi^ only 
at p ^ 1. The quasiclassical approximation^^ is valid 
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only at p\{T) ^ 1, where A(T) is the average de Brogile 
wavelength of thermally excited magnons expressed in 
lattice units. Despite the parameter e was not yet uti- 
lized for establishing of the perturbative approach valid- 
ity we suppose that this energetic quantity also is very 
important. 

In this light it seems reasonable that besides numer- 
ical estimationsiiii^ it is necessary to have for p and e 
reliable analytic formulas. The latter are just obtained 
in the present paper in the second order of perturba- 
tion theory near the exact rung-dimerized ground state 
(1). As it was already mentioned technically these cal- 
culations are possible because the first order correction 
term in the perturbative expansion for |0) contains only 
(previously obtained by the author^ilSi) two-rungon exci- 
tations around |0)r-d- At p ^ 1 and e <C 1 the obtained 
formulas agree with the corresponding numerical data^. 



II. STRUCTURE OF THE HAMILTONIAN 

It is more convenient to represent the general spin lad- 
der Hamiltonian in the following formic 



where 



TTT 

Tjd 
rrrr 

Tjll 
TTCld 



J, 



Jnorm 4\^^ 

is given by the formulas 



Sl,n+1 • 82,11+1), 
Jli^l.n ■ Si,„+i + S2,n • 82,11+1), 
Jd{^l,n ■ S2,n+1 + S2,„ ' Si,„+i), 
Jrri^l.n ' S2,ri ) (Si_„+i • S2,„+l), 

)(S2,« • S2,„+l) 
Jddi^l.n ■ S2,n+l)(S2,ri ' Si,„+i), 
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Ji = ^(2J, 



SJrr Jll J dd 

J2 = i (4( Ji - Jd) + Jll - Jddt) , 



J3 = 
J4 = 



H — Hr-n 



JeV, 
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where 



Je = ^{'^(■^i ~ Jd) - Jll + Jdd 



(13) 



Hr-d = + Qn+l) + J2(*n ' *n+l + *« ' *„+l) 

n 

+ J5(S r, (6) 



Let us notice that one can inverse the relations (8) and 
express the spin operators 



Si, 



and 



Here 
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(14) 



'S'n — 2 (Sl,n ^ S2,n) — i[Sl,„ X S2,Ti], 

*n = ^(Sl,„ - S2,„) -|-z[Si,„ X S2,„], (8) 

Sri = Si ,j + 82,71 

Qn oS^i ^71 ' ^717 (9) 



(Si.n {i ~ 1,2) are spin-1/2 operators associated with 
n-th rung). According to 

^^n^O, ^^Jlt = Sat\0)n, (10) 

the operators ^ and ^ may be considered as (nei- 
ther Bose, nor Fermi) creation-annihilation operators for 
rung-triplets. Correspondence between (5)-(7) and the 
traditional representation 

oo 

JJ = X] Hn,n+l + ^n,n + l + ^n,n+l 

n— — oo 

+ -^n,n+l + ^n,n+l + ^n,n+l + Jnorm; (H) 



As it was already mentioned^*^ the operators (8) do not 
coincide with very similar ones introduced in Refs. 3 and 
17,18 (as well as in other papers of these authors). Indeed 
the representation (14) is similar to the ones suggested in 
Refs. 17 and 3, but in fact is not identical to any of them 
because the analogs of and treated in these Refs. 
act in extended vector spaces. That is why for example 
the "inverse" formula (8) fails for these operators. The 
operators suggested in Ref. 18 in fact act in the same 
vector space as (8) but have a slightly different form. 
According to (6) 



P = 



2Ji, 
d_ 



-(Jie). 



(15) 
(16) 



Without loss of generality one can imply J2 > because 
a substitution J2 — !■ — J2 is equivalent to permutation of 
spins in all odd (or even) rungs. 

As it readily follows from (4), (6) and (10) [Q, Hr-d] = 
0. So the Hilbert space related to Hr-d splits into a direct 
sum 



m. 



(17) 
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where each Ti."^ corresponds to m-rungon sector. 

For rather big Ji vector (1) is the (zero energy) ground 
state of Hr-d- The complete set of inequahties on Ji 
which guarantees the exact rung-dimerization is not yet 
obtained except the following onf^ 



Ji - J2 > 0. 



According to (7) and (10) 



(18) 



(19) 



So the operator V destroys the exact rung-dimerization. 
However, it seems reliable that the region Ji > ^2,3,4,5 > 
Jq should contain the rung-dimerized phase for which the 
operator JqV may be treated perturbatively and 



p < 1, e < 1. 



(20) 



The parameter Ji characterizes a chemical potential 
of an exited rung (or rungon mass) while J2 its kinetic 
energy. The couplings J3 — J5 describe a spin-dependent 
rungon-rungon interaction. Finitely the parameter Jq 
governs creation-annihilation of singlet rungon pairs. 



III. EVALUATION OF THE GOVERNING 
PARAMETERS 

According to Eqs. (3) and (19) r-d{0\V\0)r-d = 0, so 

£ = £bound + £scaU + o{Jq), (21) 

where 

|2 



Abound 
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Ji j.^ \{bound\V\0)r-d\ 



2Ji N^oa NEbound 

J2 r \{q,scatt\V\0)r-d? 



SttJi Jq Escatt{q) 



dq (22) 



(we have utilized l/NJ2g ^ l/(27r) d(2g)) are con- 
tributions from singlet, translationary invariant bound 
and scattering two-rungon states^'— 



\bound) 



\q, scatt) 



A™-"...|1)^„...|1)^„, 

rn<n 

V2 



V3(A§ - 2Aocosg+ 1)7V 
X ^ a(gf,n-m)...|l)^...|l)^„. 



(23) 



Here < g < tt is a half of relative wave number 

a(g, n) ~ sinnq — Aq sin (n — l)q, (24) 
and Ao = (J3 - 2J4 + 4J5)/(2J2). 



The corresponding energies are the following 

Ebound — 4Ji + 2J2(Ao -f- — — 



= 2E. 



E, 



tiq) 



(Aq + 1)- 
2Ao 

4(Ji + J2 cosq). 



gap 



AE, 



rung ■} (25) 

(26) 



where Egap = 2(Ji — J2) and AEmng = 4J2 are the 
rungon gap and energy width^" . 

As it follows from (23) the translationary invariant 
bound state exists only for 

|Ao| > 1. (27) 

(however as it follows from (25) and (26) in a real com- 
pound it should metastable at Ag > 1). 
According to (7), (10) and (23) 

Vising 



{q, scatt\V\0)r-d 



{bound\V\0)r-d = VN 



v/(A^-2Ao cosg+1)' 
V3(Ag-l) 
Ao ■ 



(28) 



Substituting (28) into (22) and utilizing (25)-(27) one 
can readily obtain 



£ scatt 



3 P 



8JiAc 



Abound — 9(Ao — 1) 



'^Jl^Q^bound 



(29) 



{Q{x) = 1 for X > and 6(a;) = for x < 0), where 
If [z^- Ifdz 



2ni 



1^1=1 z{z - Ao)(z - l/Ao)(J222 + 2Jiz + J2) 



(30) 

For calculation of T we note that according to (18) 
the pole z = — (Ji -I- \/ — Ji)/^^ li^s outside the 
unit circle. Therefore only residues in the poles 2 = 0, 
z = {\/ Ji — J| — Ji)/ J2 and one of the appropriate 
points 1/Ao (for |Ao| > 1) or Ao (for |Ao| < 1) give 
contributions to X. Performing the calculations one can 
obtain 



3J| 



8J1J2A0 



1 



_ J2|Ag-l|+2Aov/jf~7| 
[2AoJi + (A2 + 1) J2] 

Then according to (16) 

P — Pbound + Pscatt + o{Ji), 



(31) 
(32) 



where 



Pscatt 



4[2Ao Ji + (A^ + 1) J2]^ 



(A^ + 1) Ji + 2Ao J2 



•^2 



|A^-1| 



Pbo 



6e(A2~i)($|-^ 



^oEbound 



(33) 
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The denominators a/ — J| oc ^Ega%, and E'^^^^^^ in 
(33) describe the ground state (1) destruction caused by 
condensation of scattering and bound singlet pairs. 

IV. COMPARISON WITH NUMERICAL 
RESULTS 

Here we compare the DMRG data-'^ for the case 

< Jr < 2, Jl^l, Jrr = Jll = Jdd = 0, (34) 

(also in Ref. 11 was taken Jnorm — 0, so we had to shift 
the data for varepsilon) with the corresponding results 
obtained analytically with the use of the perturbative for- 
mulas (21), (29), (31) and (32), (33) It is more instructive 
to do it in the case of maximal dimerization (considered 
in Ref. 11) namely for Jr = 2. 

For Jd = the Ref. 11 gives p = 0.11, e = 0.11 (our 
results: p = 0.22, e = 0.14). For Jd ^ 0.2, p = 0.10 
and £ = 0.08 (our results: p = 0.14, e = 0.09). For 
= 0.4, p = 0.07 and e = 0.05 (our results: p = 0.09, 
e = 0.05). For Jd = 0.6, p = 0.04 and e = 0.03 (our 
results: p = 0.04, e ^ 0.03). For Jd = 0.8, p = 0.01 and 
e < 0.01 (our results: p = 0.01, e = 0.007). Comparing 
numerical and analytic results one may conclude that the 
suggested perturbative approach works better for frus- 
trated ladders. Really at Jd = 0.8 even for Jr = 1.6 the 
numerical^i result p = 0.03 coincides with the theoret- 



ical one. Also for Jd = 0.8 and J,. = 1.5 in the both 
approaches e = 0.02. 

A discrepancy between the numerical data and analyt- 
ical predictions for non frustrated spin ladders indicates 
that besides dynamics of single rungons their properties 
are governed by creation-annihilation processes. However 
the latter may be suppressed by strong antiferromagnetic 
rung coupling. 



V. SUMMARY 

In this paper the perturbative formulas for the gas pa- 
rameter (Eqs. (32), (33)) and ground state energy energy 
per rung (Eqs. (21), (29), (31)) were obtained for spin 
ladders belonging to the rung-dimerized phase. At strong 
diagonal frustration as well as at rather strong antiferro- 
magnetic rung coupling the result agrees with the pre- 
vious DMRG calculations. In this region the conditions 
(20) are satisfied so one may conclude that the system 
really lies in the rung-dimerized phase. 

It will be interesting to compare the analytical and 
numerical approaches more precisely and in a wider range 
of the coupling parameters including four-spin terms. 

The author is grateful to M. I. Vyazovsky, S. V. 
Maleyev, S. L. Ginzburg and A. V. Syromyatnikov for 
helpful discussions. 
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